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摘要 : 分 析 现 有 的 脉冲 星 搜索 方法 对 于 即将 开展 的 FAST 望远镜 早期 脉冲 星 搜索 的 可 行 
性 和 必要 性 ， 尝 试 脉冲 星 搜索 数据 处 理 的 软 、 硬 件 加 速 。 回 顾 并 总 结 了 现 有 的 脉冲 星 搜索 方 
法 和 这 些 方法 对 应 的 发 现 ， 分 析 这 些 方 法 的 特点 。 编 写 了 基于 PRESTO 的 并 行 数据 处 理 程 
序 。 估 计 了 使 用 显卡 计算 加 速 来 完成 单 脉冲 搜索 的 可 行 性 和 效率 。 编 写 的 基于 PRESTO 的 并 
行 数 据 处 理 程序 使 得 处 理 帕 克 斯 望远镜 多 波束 巡天 的 一 个 典型 大 小 (100 MB) 的 文件 的 时 间 
从 约 95 min 缩短 到 短 于 10 min。 使 用 显卡 加 速 的 单 脉冲 搜索 实现 了 大 约 20 倍 的 加 速 ， 数 据 
处 理 时 间 短 于 获取 观测 数据 的 时 间 。 对 于 使 用 望远镜 早期 观测 数据 进行 脉冲 星 搜索 而 言 ， 用 
jm 于 短 轨 道 周 期 脉冲 星 搜索 的 相位 调制 搜索 算法 是 没有 必要 的 ， 用 于 搜索 较 长 轨道 周期 的 加 速 
搜索 算法 也 不 十 分 必要 ， 而 用 于 长 自转 周期 脉冲 星 搜索 的 快速 折 司 算法 仅 在 需要 搜索 长 自转 
周期 的 弱 脉 冲 星 时 才 需 要 使 用 。 在 数据 处 理 过 程 中 ， 软 、 硬 件 等 方面 可 以 加 快 数据 处 理 的 速 
ais 度 ， 可 以 使 数据 处 理 用 时 短 于 获取 数据 的 时 间 。 
c 关键 词 : 脉冲 星 ; 射电 望远镜 ; FAST 
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~N 到 目前 为 止 ， 已 经 发 现 超过 2500 Wika RET, ， 其 中 ， 一 些 特殊 脉冲 星 的 发 现 带 来 了 新 的 认识 。 
脉冲 双星 … 和 双 脉 冲 星 汪 的 发 现 提供 了 在 更 高 精度 上 检验 广义 相对 论 的 机 会 毫秒 脉冲 星 ” 的 发 现 
进一步 限制 了 致密 星 的 物 态 方程 ;球状 星团 中 大 量 训 秒 脉 冲 星 (例如 ， 文 [4] 发 现 Terzan 5 球状 星团 
> 中 的 21 颗 毫 秒 脉冲 星 ) 的 发 现 对 球状 星团 的 演化 过 程 给 出 了 限制 ; 银 心 附近 脉冲 星 的 发 现 表 明 ， 银 
v 心 附 近 具 有 非常 高 的 电子 密度 (如 位 于 银 心 附近 的 J1745-2900， 色 散 高 达 1778 em? pe 7) ; 在 M31 
f 视线 方向 的 X S34 cop S 的 发 现 意味 着 河 外 脉冲 星 的 科学 即将 开始 。 

更 高 效 地 发 现 新 的 脉冲 星 需要 在 硬件 或 软件 上 有 新 的 突破 。 在 硬件 方面 ， 更 大 的 望远镜 接收 面 
让 、 更 低 的 接收 机 噪声 都 有 助 于 提高 观测 的 灵敏 度 ; 而 更 强大 的 计算 机 则 可 以 提高 数据 处 理 效率 。 在 
软件 方面 ， 为 了 更 快 地 处 理 观 测 数据 ， 更 有 效 地 移 除 观测 数据 中 的 干扰 以 及 搜索 观测 数据 中 可 能 存在 
的 脉冲 星 ， 尝 试 各 种 方法 用 于 搜索 观测 数据 中 的 周期 信号 。 作 为 目前 在 建 的 全 世界 最 大 的 单 天 线 望 远 
镜 ，500 m 口径 球面 射电 望远镜 (Five-hundred-meter Aperture Spherical radio Telescope, FAST ^ ) 已 于 
2016 年 9 月 25 日 竣工 。 本 文 回 顾 过 去 十 几 年 常用 的 低频 射电 脉冲 星 观测 数据 的 处 理 方法 ， 分 析 这 些 
方法 对 于 500 m 口径 球面 射电 望远镜 早期 脉冲 星 搜索 的 必要 性 和 可 行 性 ， 并 尝试 使 用 简单 的 程序 完成 
并 行 加 速 脉冲 星 搜 索 的 处 理 流程 ， 还 尝试 了 对 搜索 过 程 进 行 图形 处 理 器 加 速 。 


1 脉冲 星 搜索 方法 


最 早 的 脉冲 星 可 以 说 是 Bell 女士 用 肉眼 发 现 的 ”， 而 数 十 年 后 的 今天 ， 脉 冲 星 的 发 现 已 经 离 不 
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开 电 子 计算 机 。 从 最 初 在 时 域 信号 中 搜索 ， 到 通过 消 色 散 和 傅 里 叶 变 换 获 得 频 域 信息 ， 再 到 今天 ， 脉 
冲 星 搜索 方法 已 经 发 展 得 相对 完善 。 
如 果 尝 试 在 脉冲 星 搜 索 中 应 用 新 技术 ， 那 么 无 论 是 使 用 原理 稍 有 改进 的 程序 ， 还 是 使 用 和 前 人 工 
作 完 全 不 同 的 周期 信号 搜索 技术 ， 通 常 的 做 法 是 基于 历史 观测 数据 集 ， 结 合 新 技术 ， 尽 可 能 多 地 找到 数 
据 集中 的 已 知 脉冲 星 ， 并 尝试 搜索 数据 集中 可 能 存在 的 、 在 之 前 的 数据 处 理 中 未 被 发 现 的 脉冲 星 。 
在 新 技术 的 尝试 过 程 中 ， 倾 向 于 使 用 帕克 斯 望远镜 多 波束 脉冲 星 巡 天 (Parkes Multibeam Pulsar 
Survey, PMPS) 数 据 作 为 基础 。 从 2000 年 的 第 1 篇 论文 1( 按 发 表 时 间 ) 到 现在 ， 共 有 25 篇 论文 记述 
了 基于 帕克 斯 望远镜 多 波束 脉冲 星 巡 天 数据 发 现 的 脉冲 星 ， 共 计 833 颗 @ 。 
根据 文 [10] 的 描述 ， 帕 克 斯 望远镜 多 波束 脉冲 星 巡 天 是 指 用 帕克 斯 望远镜 配合 13 波束 接收 机 ， 
对 银 经 260" 到 50”、 银 纬 -5° 到 + 范围 内 的 天 区 进行 的 巡天 。 此 巡天 的 第 1 篇 综述 文章 ( 文 [11]) 
发 表 于 2001 年 ， 由 于 在 观测 的 同时 已 经 开始 处 理 观 测 数据 ， 一 些 非常 重要 的 结果 已 在 此 之 前 发 表 ， 
EHR J1119-6127 FI J1814-1744 具有 很 强 的 磁场 5” ; 发 现 1141-6545 处 于 一 个 轨道 周期 只 有 5 
小 时 的 双星 系统 [1 ; 发现 1811-1736 所 处 的 双星 轨道 的 偏心 率 高 达 0.828 02 (3); URAM 
J1740-3052 的 伴星 质量 可 以 达到 11 倍 太阳 质量 。 文 [14] 指 出 ， 帕 克 斯 望远镜 多 波束 脉冲 星 巡 天 
数据 中 已 经 发 现 了 很 多 新 脉冲 星 ， 但 仍 有 不 足 ， 例 如 ， 发 现 了 非常 少量 的 轨道 周期 短 于 一 天 的 毫秒 肪 
> 冲 双星 。 文 [15] 以 一 些 同样 使 用 帕克 斯 望远镜 多 波束 脉冲 星 巡 天 数据 发 现 的 新 脉冲 星 回答 了 这 个 问 
题 。 在 文中 报告 了 新 发 现 的 11 颗 脉 冲 双 星 、15 颗 毫 秒 脉 冲 星 和 1 颗 自 转 周期 长 达 7.7s 的 脉冲 星 ; 


Co 而 用 于 搜索 的 方法 包括 日 后 被 广泛 用 于 双星 搜索 的 加 速 搜 索 方 法 (Acceleration Search) 。 到 2015 Æ, 
T 帕克 斯 望远镜 多 波束 脉冲 星 巡 天 数据 系列 论文 已 经 发 表 了 7 篇 (最 近 的 一 篇 是 文 [16] ) 。 


一 从 发 现 旋转 射电 暂 现 源 !71(Rotating Radio Transients, RRTAs) 开始 ， 大 家 在 帕克 斯 望远镜 多 波束 
Lm 脉冲 星 巡 天 数据 中 尝试 各 种 新 的 方法 ， 成 功 发 现 了 新 脉冲 星 的 方法 都 已 发 表 。 这 些 论文 总 结 在 表 1 
~N 中 。 在 这 10 篇 论文 中 ，3 篇 是 关于 旋转 射电 暂 现 源 的 搜索 ，3 篇 是 关于 搜索 结果 ( 疑似 信号 ) BG i 
; 选 排序 方法 ，1 篇 采用 了 新 的 消除 干扰 和 噪声 的 方法 ，2 篇 是 关于 全 新 的 搜索 方法 ，1 篇 采用 了 原理 
稍 有 改进 的 另 一 种 软件 ， 共 找到 了 107 颗 新 脉冲 星 。 


= 表 1 2006 年 以 来 使 用 各 种 方法 在 帕克 斯 望远镜 多 波束 脉冲 星 巡 天 数据 中 找到 新 脉冲 星 的 总 结 
PP Table 1 A summary of pulsars discovered in PMPS data since 2006 
C Paper Instruction New Pulsar Numbers 
c McLaughlin et al. 2006: Finding a new type of object, RRATs 11 
Keith et al. 2009! Using new method for candidates ranking 28 
Eatough et al. 20091”? Using new method for RFIs ( Radio Frequency Interferences , 4 


RFIs) removing 


Keith et al. 2009?! Using new method for candidates ranking 1 
Eautough et al. 2010?" Using neural network for candidates ranking 1 
Keane et al. 20107 RRATs search 10 
Keane et al. 20117? RRATs search 7 
Mickaliger et al. 2012 Using SigProc for new pulsar search 5 
Knispel et al. 20137: The project of Einstein@ Home, systematical reprocessed 24 


the data, volunteers joined for candidates identification 


Eatough et al. 201379 Coherent acceleration search 16 


@ http://www.atnf.csiro.au/ people/ pulsar/ psrcat/psrcat. help.html 
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2 不 同 的 脉冲 星 搜索 方法 对 FAST 早期 脉 促 星 搜索 的 必要 性 和 可 行 性 


脉冲 星 的 搜索 过 程 大 致 可 以 分 为 检查 并 剔除 数据 中 的 射频 干扰 (Radio Frequency Interference, RFI) , 
数据 处 理 和 搜索 (包括 如 消 色散 、 傅 里 叶 变 换 、 搜 索 周 期 信号 等 过 程 ) 以 及 疑似 信号 筛选 3 个 步骤。 

将 来 500 m 口径 球面 射电 望远镜 首先 在 低频 进行 脉冲 星 搜索 ， 观 测 频率 大 约 为 280 MHz 到 1 620 
MHz。 在 这 个 频率 范围 内 ,包含 了 电视 信和 号、 公共 通信 频率 (433 MHz) 、 手 机 通信 频率 (900 MHz) MA 
及 全 球 定位 系统 信号 (1 227 MHz, 1575 MHz) 等 干扰 信号 。 这 些 较 强 的 干扰 信号 在 消 色 散 的 过 程 中 
县 加 到 观测 数据 中 ， 对 后 续 的 数据 分 析 造 成 很 大 影响 ， 甚 至 对 于 疑似 信号 的 判断 有 本 质 的 影响 。 图 1 
是 使 用 PRESTO9 ( PulsaR. Exploration and Search TOolkit, PRESTO) 中 prepfold 命令 处 理 一 个 包含 了 已 
知 脉冲 星 信 号 的 帕克 斯 望远镜 多 波束 脉冲 星 巡 天 数据 文件 的 结果 。 一 些 通道 由 于 受 射 频 干扰 的 影响 已 
经 被 程序 去 除 (图 1 中 间 一 列 从 上 到 下 第 1 幅 灰 度 图 中 的 白色 区 域 ) ， 但 仍 有 干扰 (图 1 中 左下 灰 度 图 
中 的 斜 条 纹 ) 存 在 于 消 色 散 后 的 时 域 信 号 中 。 这 颗 脉 冲 星 信号 较 强 ， 所 以 数据 处 理 结果 没有 受到 干扰 
信号 的 显著 影响 。 帕 克 斯 望远镜 处 于 非常 宁静 的 射电 环境 中 ,仍然 面临 射频 干扰 的 问题 ， 对 于 500 m 
口径 球面 射电 望远镜 ， 虽 然 在 规划 和 建设 中 试图 避免 大 量 可 能 的 干扰 信号 ， 但 射频 干扰 的 识别 和 剔除 
在 数据 处 理 中 仍然 是 必须 的 。 


Search Information 
2 Pulses of Best Profile Candidate: ACCEL_Cand_1 RApooo= 18:34:03.0631 DECpooo7 -07:30:13.9300 
Telescope: Parkes Folding Parameters 
Epoch topo= 51 373.496 382 87037 DOF,g— 60.42. x2,4—48,621 _ P(Noise)-0 (52.50) 
Epochpary= 51 373.502 632 229 65 Dispersion Measure (DM; pc/cm?) = 287.116 
Tsample= 9.000 25 Ptopo(ms) = 512.987 6(14) Ppary(ms) = 512.977 0(14) 
Data Folded = 734 003 2 P' topo (s/s) = 0.1(5.9) X 10° P'bary (s/s) = 0.0(5.9) X 10? 
Data Avg = 49.26 P'topo (8/87) = 0.0(2.1) X 107 P"bary (s/s?) = 0.0(2.1) X 107! 
Data StdDev = 4.707 Binary Parameters 
Profile Bins = 64 Pop(s) = N/A e=N/A 
Profile Avg = 5.65e+06 a,sin(i)/c(s) = N/A o (rad) = N/A 
Profile StdDev = 159 4 
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图 1 PRESTO 处 理 包含 已 知 脉冲 星 信号 的 数据 文件 的 结果 。 数 据 中 包含 的 射频 干扰 在 图 像 中 清晰 可 见 ， 不 过 
由 于 影响 最 大 的 射频 干扰 信号 已 经 被 移 除 且 脉冲 星 信 号 较 强 ， 并 没有 对 处 理 结果 造成 本 质 影 响 
Fig.1 An example of RFI in observation data from PRESTO output. The RFI (diagonal stripe in the grey scale image at the 
lower left corner) can be seen clearly in the time domain data, but since the pulsar signal is stronger than this RFI and 


the most influential RFI is removed, this RFI does not exert major effect on the data reduction result 


@ http://www.cv.nrao.edu/ ~ sransom/ presto/ 
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搜索 脉冲 星 的 周期 信号 的 方法 主要 有 3 种 : 相位 调制 搜索 (Phase Modulation Search'*”! ) 、 加 速 
搜索 ( Acceleration Search! 5?" ) 和 快速 折 肢 算法 (Fast Folding Algorithm ^" ). ZA 500 m 口径 球面 射电 
望远镜 早期 脉冲 星 搜索 的 特点 ， 逐 个 分 析 这 些 搜索 方法 的 必要 性 和 可 行 性 。 

500 m 口径 球面 射电 望远镜 早期 观测 采用 漂移 扫描 ( Drift Scan) 的 模式 ， 将 望远镜 指向 观测 目标 即 
将 经 过 的 位 置 ， 等 待 源 从 波束 中 经 过 ， 并 在 源 经 过 波束 的 时 间 段 记录 观测 数据 。 因 此 ， 单 次 的 观测 时 
间 与 波束 大 小 有 关 ， 而 波束 大 小 与 观测 频率 有 关 。 对 于 工 波 段 ， 以 1.4 GHz 为 例 ， 波 束 宽 度 约 为 3'， 
而 源 从 波束 中 经 过 的 时 间 ， 也 就 是 单 次 的 观测 时 间 ， 大 约 为 12 s; 而 如 果 将 观测 频率 降低 到 300 MHz, 
则 波束 宽度 增加 到 14'， 观 测 时 间 增 加 到 56 s。 在 早期 观测 数据 的 处 理 中 ， 主 要 使 用 PRESTO 软件 。 

相位 调制 搜索 主要 针对 轨道 周期 非常 短 的 脉冲 双星 ， 在 观测 时 间接 近 或 长 于 双星 轨道 周期 的 情况 
下 ， 这 种 方法 可 以 获得 较 好 的 效果 。 eo dn a ne 500 m 口径 球 
面 射 电 望 远 镜 早期 观测 的 单 次 观测 时 间 非 常 短 ， 远 远 短 于 轨道 周期 。 在 望远镜 早期 脉冲 星 搜索 数据 处 
理 中 没有 必要 使 用 相位 调制 搜索 方法 。 

加 速 搜 索 采用 一 个 定常 加 速度 近似 脉冲 双星 轨道 运动 的 加 速度 ， 从 而 消除 观测 数据 中 由 于 双星 的 
轨道 运动 导致 的 脉冲 到 达 时 间 的 变化 。 加 速 搜 索 适 用 于 双星 轨道 周期 长 于 观测 时 间 ， 尤 其 是 双星 轨道 
周期 远 长 于 观测 时 间 的 情况 。 可 以 定性 估计 加 速 搜索 对 于 圆 轨道 的 脉冲 双星 是 否 有 效 : 

A. max P = 4m’xc/P,”, 
其 中 ,x 是 天 球 投 影 上 的 轨道 半 长 轴 ， 以 光 秒 为 单位 ; c 是 光速 ; P, 是 轨道 周期 ， 以 秒 为 单位 。 

A. survey. P = zmaxcP/( T^h) , 

其 中 ，zmax 是 PRESTO 中 的 加 速 搜索 设置 的 参数 ， 用 来 控制 加 速 搜索 中 尝试 的 加 速度 的 上 限 ; P 是 
脉冲 星 的 自转 周期 ; 了 是 观测 时 间 ; h 是 搜索 中 谷 加 的 谐 波 数 。 如 果 

A_max_P > A_survey_P , 
那 就 说 明 设 置 的 加 速度 (zzmox 值 ) fi], Suis ST T 18 BUR "TE SX AER E AREAS PE 47 中 
的 已 知 脉冲 双星 9， 取 望远镜 观测 时 间 为 1 min, zmax 数值 取 10, 计算 结 果 见 表 2。 


表 2 使 用 球状 星团 杜 鹏 座 47 中 的 脉冲 双星 参数 估算 500 m 口径 球面 射电 望远镜 脉冲 星 漂移 扫描 巡天 的 数据 处 理 中 加 
速 搜索 的 有 效 性 。 望 远 镜 无 法 观测 球状 星团 杜鹃 座 47， 仅 使 用 脉冲 星 的 参数 估计 搜索 的 有 效 性 
Table 2 An Estimation of FAST pulsar search efficiency of the known pulsars in the globular cluster 47 Tucanae. Because 


FAST can’t observe this globular cluster, we only use the orbit parameters of these pulsars for estimation 


PR Spin Period Orbit Period x A_survey_P A max P 
/ms day (light second) 1 minute obs. zmax 10 
J0024-7205E 3. 536 33 2. 256 84 1. 981 842 38 184. 18 0. 62 
J0024-7204H 3. 21034 2.3577 2. 152 813 3 167. 21 0. 61 
J0024-72041 3. 484 99 0. 229 79 0. 038 447 2 181. 51 1. 16 
J0023- 7203J 2. 100 63 0. 120 66 0. 040 402 22 109. 41 4. 40 
J0024- 72040 2. 64334 0. 135 97 0. 045 1512 137. 67 3. 87 
J0024- 7204P 3.64302 0. 147 2 0. 038 189. 74 2. 78 
J0024-7204Q 4. 033 18 1. 189 08 1. 462 209 4 210. 06 1. 64 
J0024-7204R 3. 480 46 0. 066 2 0.0334 181. 27 12. 09 
J0024— 72048 2. 83041 1. 201 72 0. 766 270 3 147. 42 0. 84 
J0024-7204T 7. 588 48 1. 126 18 1.338 502 395. 23 1. 67 
J0024-7203U 4. 342 83 0. 429 11 0. 526 950 1 226. 19 4. 54 
J0024-7204V 4.81 0. 227 16 0. 8513 250. 52 26. 17 
J0024-7204W 2.35234 0. 133 0. 243 5 122. 52 21. 84 
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在 1 min 的 观测 数据 中 ， 这 些 脉冲 星 的 轨道 运动 导致 的 脉冲 到 达 时 间 对 于 搜索 的 影响 很 小 。 一 方 
面 ， 这 些 脉 冲 星 都 是 毫秒 脉冲 星 ， 自 转 非 常 快 ， 对 于 一 般 的 脉冲 星 ， 比 如 自转 周期 为 0. 1s, A_survey_P 
的 数值 会 进一步 增 大 ， 这 意味 着 可 以 使 用 更 小 的 zmax 值 ， 比 如 仅 取 1 或 者 不 进行 加 速 搜 索 (zmax=0， 
则 相当 于 不 进行 双星 搜索 ， 但 是 仍然 有 可 能 找到 脉冲 双星 的 信号) 。 

对 于 PRESTO 的 加 速 搜索 程序 来 说 ，10 甚至 更 低 的 zmax 数值 ， 消 耗 的 计算 时 间 和 内 存 都 非常 
少 。 因 此 ， 对 于 500 m 口径 球面 射电 望远镜 早期 脉冲 星 搜索 来 说 ， 可 以 不 采用 加 速 搜 索 ， 如 果 需 要 搜 
索 轨 道 周期 特别 短 的 脉冲 双星 ， 可 以 尝试 加 速度 较 小 的 加 速 搜索 。 

快速 折 笃 算法 主要 用 于 搜索 自转 周期 特别 长 的 脉冲 星 。 尝 试用 PRESTO 中 的 prepfold 处 理 包 含 了 
自转 周期 长 达 7.7s 的 脉冲 星 1001-5939 7" 的 帕克 斯 望远镜 多 波束 脉冲 星 巡 天 数据 文件 ， 其 信号 的 
信 噪 比 可 以 达到 69. 3( PM0081_02471.sf) ， 但 使 用 PRESTO 搜索 只 能 得 到 这 个 信号 的 二 次 谐 波 ， 信 噪 
比 为 51. 8。 显 然 ， 如 果 在 搜索 中 得 到 的 并 非 是 脉冲 星 信 号 的 基 频 信号 ， 而 是 谐 波 信号 ， 那 么 搜索 出 
的 谐 波 信号 的 信 噪 比较 基 频 信号 有 所 降低 。 在 这 种 情况 下 ， 较 强 的 长 周期 脉冲 星 仍 然 可 能 被 搜索 到 ， 
但 是 也 会 导致 一 些 弱 信 号 被 遗漏 。 因 此 ， 如 果 要 搜索 长 自转 周期 的 弱 脉 冲 星 ， 有 必要 尝试 快速 折 秋 算 
法 处 理 观测 数据 。 

PRESTO 中 不 包含 可 以 进行 快速 折 秋 算法 的 程序 ， 要 使 用 这 一 算法 ， 需 要 使 用 Sigproc® 中 的 ffa 
程序 。 

另外 , 文 [33] 采 用 分 段 搜索 (Segmented Search) 的 方法 累加 了 对 于 球状 星团 杜 鹏 座 47 的 大 约 
1 100 小 时 的 观测 数据 并 成 功 找到 了 新 的 脉冲 星 。 对 于 500 m 口径 球面 射电 望远镜 早期 观测 中 单 次 观 
测 每 个 源 的 时 间 有 限 这 一 情况 ， 采 用 分 段 搜索 的 方法 ， 受 加 不 同时 间 段 的 观测 数据 ， 提 高 探测 灵敏 度 
是 可 行 的 。 这 一 方法 可 以 用 于 今后 500 m 口径 球面 射电 望远镜 脉冲 星 搜索 中 。 需 要 指出 ， 文 [33] 搜 
索 的 是 球状 星团 中 的 脉冲 星 ， 所 以 仅 使 用 了 一 个 色散 值 ， 而 对 于 一 般 的 搜索 ， 需 要 尝试 数 百 个 色散 
值 ， 因 而 计算 量 也 增加 数 百倍 。 

在 疑似 信号 的 筛选 阶段 ， 显 然 使 用 计算 机 程序 ， 无 论 是 简单 的 筛选 ， 比 如 按照 信 噪 比 排序 ， 还 是 
使 用 复杂 的 方法 ， 比 如 人 工 智能 ， 都 是 可 行 的 。 但 这 些 程序 或 许 会 遗漏 非常 少量 的 真实 的 脉冲 星 信 
号 。 也 可 以 结合 人 工 ， 比 如 Einstein@ Home? 。500 m 口径 球面 射电 望远镜 每 年 可 能 产生 数 以 亿 计 的 
脉冲 星 疑似 信号 。 结 合计 算 机 和 人 工 可 能 是 更 好 的 选择 。 对 于 海量 的 疑似 信号 ， 首 先 使 用 计算 机 程 
序 ， 将 其 中 质量 不 好 的 疑似 信号 剔除 ， 再 人 工 对 剩余 的 疑似 信号 作 判 断 。 


3 PRESTO 脉冲 星 搜索 加 速 和 Heimdall 单 脉 冲 搜索 加 速 的 尝试 


PRESTO 的 脉冲 星 搜索 分 步 进行 。 按 照 操作 的 先后 顺序 ，PRESTO 中 的 不 同 程序 ， 先 后 完成 分 析 
数据 中 的 射频 干扰 并 标记 ， 逐 个 消 色 散 并 输出 文件 ， 将 消 色 散 后 的 文件 逐个 进行 傅 里 叶 变换 得 到 频 域 
言 息 ， 逐 个 搜索 不 同色 散 值 的 频 域 数 据 ， 疑 似 信号 排序 以 及 生成 最 终 供 判 断 信 号 质量 的 图 表 。 

考虑 到 对 于 不 同 的 色散 数值 ， 消 色散 、 傅 里 叶 变换 、 搜 索 以 及 最 终生 成 图 表 可 以 同时 进行 ， 自 主 
编写 程序 将 这 些 过 程 并 行 化 ， 并 行 化 的 加 速效 果 非 常 显著 。 在 一 台 配 备 了 双 路 X5690( 共 12 核 24 HK 
程 ) 和 96 GB 内 存 的 计算 节点 上 的 测试 结果 如 图 2。 对 于 单个 帕克 斯 望远镜 多 波束 脉冲 星 巡 天 数据 文 
件 ，12 个 线程 同时 处 理 ， 用 时 已 经 下 降 到 大 约 1070s, 使 用 大 约 100 甚至 更 多 线程 时 ， 处 理 时 间 降 
到 900 s 以 下 。 如 果 不 采用 并 行 化 处 理 ， 单 个 数据 文件 的 处 理 时 间 超 过 90 min。 并 行 化 处 理 节 约 了 超 
过 80% 的 时 间 。 在 配置 更 高 的 计算 节点 上 ， 处 理 时 间 可 以 缩短 到 600 s 以 内 。 

未 来 500 m 口径 球面 射电 望远镜 脉冲 星 观测 数据 也 将 在 配备 了 多 核 处 理 器 的 大 型 计算 集群 上 
成 。 在 最 初 阶段 ， 将 一 些 数 据 处 理 过 程 简单 并 行 化 ， 非 常 有 利于 提高 处 理 效率 。 


at 


®©  https://sourceforge.net/projects/sigproc/ 


fo Le XC É DN. i 
^! ^V 1 4 
C hinaA IVA 


1 期 WIRT. 脉冲 星 搜索 技术 及 FAST 望远镜 脉冲 星 搜索 展望 13 


另外 尝试 单 脉冲 搜索 的 加 速 。 单 脉冲 是 指 带 有 色散 的 单个 较 强 的 脉冲 ， 它 可 能 来 自 于 已 知 的 脉冲 
星 ， 旋 转 射电 暂 现 源 ， 也 可 能 是 快速 射电 暴 。 如 果 采 用 PRESTO 进行 搜索 ， 其 步 又 和 脉冲 星 搜 索 一 
样 ， 只 不 过 在 搜索 这 一 步 中 使 用 另 一 个 专门 用 来 搜索 此 类 信和 号 的 python 程序 。 而 作为 对 照 ， 使 用 
heimdall©'*’ 。Heimdall 是 基于 NVIDIA 硬件 和 CUDA 的 图 形 处 理 器 ( Graphics Processing Unit, GPU) 加 
速 的 程序 。 经 过 测试 ， 在 一 台 安 装 了 NVIDIA M2090 计算 卡 的 服务 器 上 ， 对 于 帕克 斯 望远镜 多 波束 脉 
冲 星 巡 天 数据 ，Heimdall 只 需要 32-34 s 即 可 完成 搜索 ， 且 能 搜索 出 其 中 来 自 于 旋转 射电 暂 现 源 的 单 
脉冲 信号 。 如 果 PRESTO 搜索 需要 600s, Heimdall 使 用 图 形 处 理 器 实现 了 约 20 倍 的 加 速 。 
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图 2 使 用 不 同 线程 数目 时 单个 文件 的 处 理 时 间 。 横 轴 是 最 大 并 行 线程 数目 ， 纵 轴 是 总 用 时 。 所 有 测试 在 
同一 个 计算 节点 上 和 运行， 使 用 同一 个 文件 。 对 于 这 个 计算 节点 来 说 ， 当 线程 数 在 大 约 100 以 下 时 ， 
增加 线程 数 对 于 缩短 处 理 时 间 有 显著 效果 

Fig.2 Data reduction time decreases with increasing thread number. The X-axis is the maximal number of thread 
running at the same time; the Y-axis is the total data reduction time (in the unit of second). In this case, 


when the thread number is below 100, increasing the threads is efficient to reduce the data reduction time 


帕克 斯 望远镜 多 波束 脉冲 星 巡 天 数据 的 单个 数据 文件 的 时 间 分 辩 率 是 250 ps, FE 96 个 通道 ， 观 
测 时 长 约 为 2 100s。 而 500 m 口径 球面 射电 望远镜 如 果 采 用 100 hs 的 时 间 分 辨 率 ，!1 024 个 通道 ， 观 
测 时 长 1 min。 如 果 也 采样 1 bit， 那 么 望远镜 单 次 观测 的 数据 量 大 约 是 帕克 斯 望远镜 多 波束 脉冲 星 巡 
天 数据 单个 文件 的 三 分 之 一 。 因 此 ， 在 单 计算 节点 上 用 PRESTO 完成 脉冲 星 搜 索 的 时 间 大 约 是 200 s。 
只 要 配备 超过 4 个 节点 ， 就 可 以 实现 处 理 时 间 短 于 数据 获取 时 间 。 如 果 使 用 Heimdall 搜索 单 脉冲 ， 则 
用 时 大 约 为 10 s， 已 经 短 于 观测 时 间 。 也 就 是 说 ， 一 个 小 规模 的 集群 可 以 实现 观测 数据 的 实时 处 理 。 


4 总 结 和 展望 


在 分 析 前 人 使 用 的 一 些 搜索 脉冲 星 的 方法 的 基础 上 得 出 结论 ,在 500m 口径 球面 射电 望远镜 早 
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期 脉冲 星 搜索 中 不 需要 使 用 相位 调制 搜索 方法 ， 只 在 搜索 自转 周期 和 轨道 周期 都 非常 短 的 脉冲 双星 时 
才 需 要 使 用 加 速 搜 索 方 法 。 这 两 种 搜索 方法 都 可 以 使 用 PRESTO 软件 实现 。 而 在 计算 资源 充足 时 ， 也 
可 以 考虑 使 用 SigProc 中 的 fa 程序 ,使 用 快速 折 秋 算法 ,试图 搜索 数据 中 可 能 存在 的 自转 周期 非常 长 
的 脉冲 星 。 

尝试 了 PRESTO 数据 处 理 的 并 行 化 ， 将 处 理 时 间 缩 短 到 原先 的 大 约 40%。 也 尝试 用 于 搜索 单 脉 
冲 的 图 形 处 理 器 程序 ， 实 现 了 大 约 20 倍 的 计算 加 速 。 这 些 搜索 计算 加 速 的 尝试 有 助 于 在 今后 提高 望 
远 镜 脉冲 星 观测 数据 的 处 理 速度 。 

配合 这 些 分 析 和 加 速 的 尝试 ， 可 以 大 致 实现 500 m 口径 球面 射电 望远镜 漂移 扫描 数据 的 处 理 速 度 
快 于 数据 获取 的 速度 。 这 些 数据 处 理 包括 脉冲 星 搜索 、 单 脉冲 搜索 和 搜索 到 的 疑似 信和 号 的 识别 。 

FB 500 m 口径 球面 射电 望远镜 非常 高 的 灵敏 度 和 脉冲 星 搜索 技术 ， 在 望远镜 早期 漂移 扫描 的 数 
据 中 找到 一 定数 量 的 新 脉冲 星 。 在 数据 处 理 过 程 中 ， 建 议 首先 快速 搜索 数据 中 的 单 脉 冲 ， 结 合 图 形 处 
理 需 加 速 ， 搜 索 可 以 在 非常 短 的 时 间 内 完成 。 搜 索 单 脉 冲 的 同时 ， 使 用 PRESTO 不 设置 加 速度 值 或 者 
只 使 用 非常 小 的 加 速度 范围 搜索 其 中 的 周期 信号 。 最 后 ， 可 以 考虑 使 用 SigProc 重新 处 理 数 据 ， 主 要 
使 用 ffa 程序 ， 即 使 用 快速 折 秋 算法 搜索 其 中 可 能 存在 的 弱 的 长 周期 脉冲 星 信 号 。 

致谢 : 感谢 美 家 射电 天 文 台 的 Scott Ransom 教授 在 PRESTO 使 用 方面 的 帮助 ， 德 国 马克 思 普 
朗 克 射电 天 文 研 究 所 的 Paulo Freire 教授 在 计算 加 速 搜索 有 效 性 方面 的 协助 以 及 加 州 大 学 伯克利 分 校 
的 Dan Werhimer 教授 在 使 用 Heimdall 方面 的 指点 。 同 时 感谢 黔 科 合作 人 才 团 队 的 帮助 和 支持 。 
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Abstract; We aim to see the feasibilities and necessities of applying current pulsar search methods to 
FAST pulsar search. We test pulsar search data processing with both software and hardware acceleration. We 
discuss methods used in current pulsar search, discoveries obtained by these methods and their characteristics. 
Based on PRESTO, we write scripts to run the pulsar search processing in a parallel manner. We evaluate the 
feasibilities and necessities of GPU acceleration in single pulse search. With the scripts, the time of processing 
a ParkesMulti-beam Pulsar Survey data file ( with a typical size of 100 MB) has been shorten from approximately 
95 minutes to less than 10 minutes. Our benchmark shows that the single pulse search in one file has been 
accelerated by 20 times. The single pulse search can be done in shorter time compared with the data aquisition 
process. Used to search for binary pulsars with short orbit periods, phase modulation search is not necessary for 
FAST pulsar search; similarly, used to search for binary pulsars with relative long orbit periods, the 
acceleration search is not important for FAST binary pulsar search and the fast folding algorithm is only 
necessary when searching for faint pulsars with long spin periods. For data processing, both software and 
hardware acceleration can and should be used to reduce time cost. 
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